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ABSTRACT 

We show that one can measure the effects of microlens parallax for binary microlensing events with three well- 
measured peaks - two caustic crossings plus a cusp approach, and hence derive the projected Einstein radius te- 
Since the angular Einstein radius 9e is measurable from finite source effects for almost any well-observed caustic 
crossing, triple-peak events can yield the determination of the lens mass M = c~ 2 : (4G)fE0E- We note that, to a 
certain extent, rotation of the binary can mimic the effects of parallax, but it should often be possible to disentangle 
parallax from rotation by making use of the late-time light curve. 

Subject headings: gravitational lensing 



1. Introduction 

When microlensing experiments were proposed (Paczyriski 1986), it was clear that the lens masses would not be measured 
individually on an event-by-event basis because, among the routine observables, only the Einstein timescale ts is related to the 
mass, and only indirectly to it; 

tv = ^: h = J^- (1) 
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Here D IC \ = AU/7r rc i, 7r rc i is the lens-source relative trigonometric parallax, /j, le \ is the relative proper motion, 6*e is the angular 
Einstein radius, and Rsch is the Schwarzschild radius of the lens mass. Rather, it would be necessary to interpret the results 
statistically, to infer an average mass scale of the lenses from the mean timescale of the observed events together with statistical 



estimates of 7r rc i and fi rc \. Note that Z? rcl 1 = D L 1 — D s , where £>l and D$ are the distances to the lens and the source. 



Gould (1992) showed that it was possible to extract the individual lens mass M if one could measure two "observables": the 
angular Einstein radius and the projected Einstein radius ?e; 



^E = #E-Drcl — V 2i?Sch A'd , (2) 

M = ^r E 9 E . (3) 
The microlens parallax tte re-expresses in a convenient way. It is a two-dimensional vector whose magnitude is given by 

= -7- = — , (4) 

#E »'E 

and whose direction is that of the lens-source relative proper motion. Measurement of 9e requires that the Einstein ring be 
compared with some "angular ruler" on the plane of the sky, while to measure fs, one must compare the Einstein ring with some 
"physical ruler" in the observer plane. 

The only angular ruler to be used for this purpose to date is the angular radius of the source 8* which can be estimated from 
its de-reddened color and magnitude, and the empirical color/surface-brightness relation (eg. van Belle 1999). When the source 
passes very close to or directly over a caustic (zone of formally infinite magnification) in the lens magnification pattern, the finite 
size of the source affects the microlensing light curve. By analyzing this finite-source affected light curve, one can measure p» 
= 9* /6e, and thus 6>e as well. While this idea was originally proposed for point-mass lenses, which have point-like caustics 
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(Gould 1994; Nemiroff & Wickramasinghe 1994; Witt & Mao 1994), in practice it has been mainly used for binary lenses, which 
have line-like caustics and hence much larger cross sections (Alcock et al. 1997, 2000; Albrow et al. 1999a, 2000, 2001; Afonso 
et al. 2000). 

Typically te is of order several AU, and so the natural scale of the physical ruler must be ~ 1 AU. This is the case for the 
two main methods proposed to measure te- First, a satellite in solar orbit would see a substantially different event from the one 
observed on Earth, and from the difference, one can determine f-E (Refsdal 1966; Gould 1995). Second, the reflex of the Earth's 
orbital motion induces a wobbling of the source's passage through the Einstein ring, which in turn perturbs the light curve, and 
Te can be determined from this perturbation (Gould 1992). To date, this latter method is the only one by which ?e has been 
measured (Alcock et al. 1995; Mao 1999; Soszyhski et al. 2001; Bond et al. 2001; Smith, Mao, & Wozniak 2001; Mao et al. 
2001; Bennett et al. 2001). Unfortunately, this method requires that the event be rather long, t-g > 90 days. To the extent that the 
Earth's motion can be approximated as rectilinear, it is not possible even in principle to detect this effect. Over the relatively short 
timescales of typical events t E ~ 20 days, on the other hand the Earth's motion can be approximated as uniform acceleration 
which can lead to a potentially detectable asymmetry in the light curve (Gould, Mirala-Escude, & Bahcall 1994). However, this 
would yield a measurement of only the projection of tte, along the direction of the Earth's acceleration vector at the peak of the 
event. Hence, even if this asymmetry were detected, it would give only an upper limit on te- To measure both components of 7Te 
requires that the event last long enough for the Earth's acceleration vector to change substantially while the source lies within the 
Einstein ring. 

However, it is possible to make parallax measurements using much shorter baselines. During a caustic crossing of a binary 
lens, the flux can change by several magnitudes as the caustic transits the star. The source radius projected onto the observer 
plane #»Z? re i can be of order <~ 100 R© for a main-sequence star. Hence, observers on two continents could see fluxes differing 
by a few percent. This would again yield only one projection of 7Te, but by repeating this procedure for the two caustic crossings 
of a single event, or by observing one crossing from three continents, one could obtain a full measurement of parallax (Hardy & 
Walker 1995; Gould & Andronov 1999). An important feature of this approach is that it also allows to measure 9e, which can 
be obtained from any caustic crossing. Hence, if this method were ever applied, it would yield the mass through equation (3). 
Unfortunately, the method requires extraordinary cooperation from the event itself, from the observatory directors, and from the 
weather. Consequently, it has not to date been successfully carried out. 

It is expected that the use of the high-precision interferometric astrometry, most notably the Space Interferometry Mission 
(SIM) will be the best way to measure the individual lens mass (Paczyhski 1998; Boden, Shao, & van Buren 1998). With 10- 
yuas level astrometry of SIM, one can directly measure 7r re i- Furthermore, by monitoring the displacements in the light centroid 
caused by microlensing, 6e can be also directly measured (H0g, Novikov, & Polanarev 1995; Miyamoto & Yoshii 1995; Walker 
1995). Unlike most techniques mentioned earlier, neither of these measurements requires the event to be of a special class such as 
long time-scale or caustic-crossing events. That is, individual lens masses can be measured for common point-source/point-lens 
events. 

Here, we propose a new method to measure parallax for a subclass of caustic crossing binaries. The method requires the 
event to have at least three separate well-defined photometric peaks so that the source position can be precisely located relative 
to the lens geometry at three distinct times. Two of the peaks should be caused by caustic crossings: an entrance and an exit. It 
would be best for the third peak to be a cuspy-caustic crossing, but the chance of two fold-caustic crossing accompanied by an 
additional cusp crossing is low. In the following, we investigate events where the third peak is due either to a cusp crossing or to 
a generic cusp approach. 



If the event has two well-observed caustic crossings, plus a well-observed "cusp approach", these features then provide 
nine "empirical parameters" - parameters that can be measured directly from the light curve, without appealing to a model that 
depends upon a global lens geometry. They are the maximum flux F max , the half-duration At, and the time at the maximum i max 
of each of the three peaks. The uncertainties in the measurements of these quantities are approximately 



2. Overview 
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where <7 p h is the mean fractional error of each of the TV photometric measurements taken over the bump. Since At is generally 
of order a day, the errors in At and t max measurements can easily be of order a minute or smaller. The fact that six independent 
times can be measured with such precision, typically four or five orders of magnitude smaller than the characteristic timescale of 
the event t E , is what makes the parallax measurement feasible. 

In principle, any localizable feature found in the light curve may provide similar empirical parameters. What then makes 
peaks more significant as information posts than other features such as minima or points of inflection? The answer to this 
question is that the high precision measurements attainable for peaks can rarely be achieved for empirical parameters associated 
with the other localizable features, so that their merit as signature beacons is much weaker than for peaks. For instance, minima 
in microlensing light curves typically have widths that are ten times greater than peaks do (and by definition less flux as well). 
From equation (5), a(t) oc \J At/F, so that it is substantially more difficult to localize minima than peaks. Here, we assume 
Poisson noise (cr p h = V~F /F), and the same sampling frequency (i.e. N oc At). Points of inflection are even more difficult to 
localize than minima. 

Ordinarily, under the approximation that the lens-source relative motion is rectilinear, to specify a binary-lens light curve 
requires seven global (or geometric) parameters: d, the binary separation in units of 6>e; q, the binary mass ratio; a, the angle of 
the source-lens relative motion relative to the binary axis; ts, the Einstein timescale; uo, the minimum angular separation between 
the source and the binary center in units of #e; to, the time at this minimum; p*, the source size in units of 8e- (In addition, 
to transform the magnification curve to the specific photometric system of the observations, one also needs limb-darkening 
parameters for each wave band of observations; plus the source flux F s and background flux Ft,, for each telescope and wave 
band.) However, since the actual motion is not rectilinear, these seven parameters will not be adequate to describe the event for 
very high precisions, and in particular, subtle inconsistencies will be introduced among the nine precisely measured quantities 
mentioned earlier. We now show how these inconsistencies can lead to a parallax measurement. 



3. Measurement of Parallax 



In practice, the parallax will be measured by multi-dimensional fitting and subsequent x 2 minimization. However, it is 
instructive for two reasons to identify in a systematic way the features of the event that permit 7Te to be measured. First, this 
enables one to predict when an event will have a measurable 7Te- Second, there are technical difficulties associated with x 2 
minimization, and these can be ameliorated if the model parameterization is modified to reflect the underlying physics (Albrow 
et al. 1999b). 

To understand how 7Te is measured, we first show how some of the nine empirical parameters are related to one another in 
the absence of parallax, i.e., 7Te = 0. We will initially assume that the lens geometry (d, q) and the fluxes F s and Ft, are known 
a priori. Then, at the peak of the cusp approach i max ,ca, the source position within the Einstein ring it ca can be determined very 
precisely from the precise measurements of t max . ca and F max ca . That is, at this time t max . ca , the source must be somewhere along 
the cusp-approach ridge line and its position on that line is determined from the inferred magnification A = (F max ca — F\>)/F s . 

Now consider the angles at which the source crosses each caustic line, (pi and fa- The half-duration of a caustic crossing is 
given by At cc = t^p* esc <f>. If this motion is assumed to be rectilinear, then t^ becomes a well-defined quantity over the whole 
duration of the event, and thus, 

CSC (f> 2 At cc o 
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The source trajectory is then the straight line that passes through it ca and has caustic-crossing angles that satisfy equation (6). 
From that equation, this angle can be determined with a precision of 
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which can be estimated using equation (5). 

Next, we turn to the times of the two caustic crossings to.ee. l and to.ee, 2, at which the center of the source lies on the caustic 
lines. These are not the same as the times of peak flux t ma x,cc,i an d £max,cc,2, but they can be determined to the same precision 
(see eq. [5]) and are more convenient to work with. Let iz C c,i an d u cc , 2 be the positions in the Einstein ring of the two caustic 
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crossings. For rectilinear motion, these satisfy the vector equation, 

^cc,2 ^*cc,l ^ca ^cc.l 



^0,cc,2 io.cc. 1 ^max.ca ^0,« 



(= IMeI) ■ (8) 



Here, |/x E | =t% 1 - However, if the acceleration due to the Earth's orbital motion is not negligible, equation (8) will not in general 
be satisfied. This can be quantified in terms of St 2 , the difference between the measured value of i ,cc,2 and the one that would 
be predicted from the other measured parameters on the basis of equation (8), 

8t 2 = (io,cc,2 — £o,cc,l) — -jfimaxxa — £o,cc,l) ■ (9) 

l^ca ^cc.l | 

To evaluate the relation between St 2 and 7Te, we make the approximation that the Earth's acceleration vector projected on the sky 
a© t _L is constant for the duration of the caustic crossings and cusp approach. We first note that the magnitude of this acceleration 
is related to the parallax by 



2 



^^(^J ' (10) 



r-E 

where ip is the angle between the lines of sights towards the Sun and the event from the Earth r® is the distance between the Sun 
and the Earth, during the caustic crossings, and f2® = 2ir yr _1 . Then, after some algebra, we find 



St 2 



8in#E f-9®\ {-7r E ,||(t2-ti)+7rte,±[(t3-t2)cot^2-(t 3 -ti)cot0i]}, (11) 



h-t 2 2 \r @ /AUj 

where (tte,\\ , tte,_l) are the components of 7Te parallel and perpendicular to a ffi ,i, and where we have used the simplified nota- 
tions, <o,cc.i — * U and fmax.ca —* h- Hence, by measuring St 2 one can determine a particular projection of 7Te whose components 
are given by equation (11). 

However, since 7Te is a two-dimensional vector, measurements of two independent components are required for its complete 
determination. A second constraint is available from the width of the cusp approach At ca . We define 8 In Ai ca in analogy to 8t 2 
and after some more algebra we find, 

8 In At ca = JiE^g- {tte.h [(t 3 - h) cot fa + (t 2 - h) cot fa - (t 3 - h) cot fa] 

cot fa — cot fa \r©/AU/ 

+^E,±[(h - h) cot fa cot 0! - (t 3 - t 2 ) cot ^ 3 cot fa - (t 2 - ti) cot ^>2 cot 0!]} . (12) 

Measurement of d In Ai ca therefore gives another projection of 7r E . Furthermore, it is possible to obtain a third constraint on 7r E 
from the behavior of the light curve after the source has left the caustic but while it still remains within the Einstein ring. For 
rectilinear motion, the Einstein time scale is well determined (cf. eq. [8]). If the late-time light curve drops off faster or slower 
than indicated by this timescale, it implies that the Earth's acceleration vector has a component aligned with the direction of lens- 
source relative motion. The effect is similar to the one identified by Gould et al. (1994) but can be measured more easily because 
both t~E and to are determined very precisely from the light curve around the caustic-crossing region. The relative orientation 
of these three constraints depends on the details of the lens geometry. In principle, they could all be roughly parallel, but this 
is unlikely; in general, it should be possible to combine the three projections to measure both components of 7Te- If the event 
is sufficiently long (typically t E > 60 days), so that the Earth moves > 1 radian during an Einstein timescale then the late time 
behavior of the light curve will give information about both components of 7r E . In this case there would be the fourth constraint. 

The exact expression for the errors in the components of 7Te obtained from the first two constraints can be derived from 
equations (11) and (12), but these are extremely complicated and for that reason not very interesting. However, reasonable 
estimates of these errors can be made as follows. First, we note that the error from applying the 8t 2 constraint is dominated by the 
problem of determining the change in a from the parallax to the non-parallax case. The error in the measurement of a is given 
by equation (7). On the other hand, the change in a (chosen for definiteness to be the angle at the time of the cusp approach) is 
given by 

*28t 

(cot fa — cot fa)8a = . (13) 

H - t 2 

Hence the fractional error in 7TE,at 2 , the projection of 7Te measured by this constraint, is of the order of 

0-(7TE,<5t 2 ) <y{a) A^ 1/2 £Tph 
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where we have made use of equations (5) and (7). For typical bulge parameters, | sin^l ~ 0.2 and 7r E ~ 0.1, and with good 
photometric coverage of the caustic crossings, N~ 1 ^ 2 a p h <~ 10~ 3 , the fractional error is ~ 0.17 x (*e/50 days) -1 x [(t 2 — 
ti)/20 days] -1 and hence 7r E ^ t2 should plausibly be measurable. For the other constraint one finds similar expressions. 

We now relax our assumption that d, q, F s , and are known a priori. Actually, if d and q are known, F s and F^ can be 
easily determined from, for example, the baseline flux F s + F^, and the minimum flux observed inside the caustic A m [ n F s + F\>. 
Here ^4 m i n is the magnification at minimum which, for fixed (d, q), is virtually independent of the minor adjustments to the 
trajectory due to parallax. However, it is still necessary to relax the assumption that d and q are known. In fact, these must 
be determined simultaneously with the parallax because changes in (d, q) can have effects on the relative times of the caustic 
and cusp crossing and on the crossing angles, just as parallax can. Nevertheless, there are numerous other constraints on (d, q) 
coming from the overall light curve, and so while (d, q) and 7r E can be expected to be correlated, they should not be completely 
degenerate. Hence, even allowing for degradation of the signal-to-noise ratio for the 7r E measurement due to the correlations 
between the projections of 7r E and (d, q), it should be possible to measure 7r E with reasonable precision for events with two 
well-covered caustic crossings and a well-covered cusp approach. 



The measurement of the parallax discussed in the previous sections essentially relies upon the failure of the rectilinear ap- 
proximation of the lens-source relative motion when the trajectory is over-constrained by available observations. For a relatively 
short time scale, what is actually measured from this is an instantaneous acceleration on the source trajectory and this acceleration 
may contain significant contributions by other effects, a notable example of which is the binary rotation around its center of mass. 
The rotational motion of the binary lens projected onto the plane of the sky is observable in terms of a contraction or expansion 
of the binary separation d, and the lateral rotation of the binary axis with respect to the fixed direction on the sky w. The apparent 
observable result of the latter effect may heuristically be understood as a centripetal acceleration ~ u c uj 2 on the source motion 
relative to the (static) lens system, where u c is the angular extent of the caustic in units of the Einstein ring. For a face-on circular 
binary orbit observed at the ecliptic pole, the ratio of the Earth's acceleration to the projected acceleration due to the binary orbit 
is then 



Thus, depending on these parameters, either parallax or binary rotation could dominate. In addition, for the general case, both 
the Earth and binary acceleration would be reduced by possibly very different projection factors (while the projection factor for 
the binary rotation is from the orbital inclination of the binary, the parallactic projection is due to the angle between the direction 
to the ecliptic pole and the line of sight to the event). Although the effect of d is in general not expressible analytically, one can 
expect it to be of similar order of magnitude to the one of cj. 

Unambiguous determination of the parallax therefore requires the measurement of the acceleration at least at two different 
times, or not less than four independent constraints on the projection of the acceleration at different times. We argued above that 
there are generically at least three constraints for the types of events under discussion and that there is a fourth constraint for 
sufficiently long events. In this latter case, parallax can be unambiguously discriminated from rotation. But even when the event 
is short, so that there are only three constraints, it may still be feasible to measure the parallax. We first note that the rotation 
cannot significantly affect the late-time light curve - i.e., i E is not influenced by rotation - for it does not affect separations 
between the source and the binary center of masses. Hence, if one component of 7r E is measured from the late-time light curve, 
it may be possible using equation (15) and Kepler's Third Law, to show that parallax is more important than rotation, in which 
case the caustic crossings can be used to determine the full parallax 7r E . 

In brief, there is good reason to hope that events with three well-defined bumps can yield parallax measurements, and hence 
mass measurements. Whether this will be possible for any particular such event can only be determined by detailed modeling. 
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4. Effect of Binary Rotation 
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